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We study transport properties of a charge qubit coupling two chiral Luttinger liquids, realized by 
two antidots placed between the edges of an integer v — 1 or fractional v — 1/3 quantum Hall bar. 
We show that in the limit of a large capacitive coupling between the antidots, their quasiparticle 
occupancy behaves as a pseudo-spin corresponding to an orbital Kondo impurity coupled to a chiral 
Luttinger liquid, while the inter antidot tunnelling acts as an impurity magnetic held. The latter 
tends to destabilize the Kondo fixed point for the v — 1/3 fractional Hall state, producing an 
effective inter-edge tunnelling. We relate the inter-edge conductance to the susceptibility of the 
Kondo impurity and calculate it analytically in various limits for both v — 1 and v — 1/3. 


Introduction - Fractional quantum Hall (FQH) sys¬ 
tems [I] are strongly correlated topological states, real¬ 
ized in clean two-dimensional electron gases under a large 
perpendicular magnetic held, where the bulk contains an 
incompressible fluid and the low energy dynamics is con¬ 
trolled by chiral Luttinger liquids at the edges [2]. There 
has recently been a renewed interest in these systems 
due to a promise of the celebrated topological quantum 
computation using non-abelian anyons mm and also in 
connection to impurities in helical liquids of the quantum 
spin Hall systems [THEE]. However, there is still a consid¬ 
erable gap between theoretical and experimental studies 
of abelian anyons in the FQH edge states which moti¬ 
vates a more thorough study of their properties. Here, 
we study the problem of elastic co-tunnelling of Laugh- 
lin quasiparticles through two antidots and show that in 
certain limits, it maps to a Kondo impurity da nm em¬ 
bedded between two chiral Luttinger liquids [14, 15, 17 
mm and exhibits interesting transport signatures. 

Transport through antidots in the FQH regime has 
been studied in the past, experimentally [2QH251 and the¬ 
oretically [26ff29] in a regime where the transport was 
dominated by correlated but incoherent transfers of in¬ 
dividual quasiparticles. In contrast, in this paper we are 
interested in a regime where this sequential tunnelling is 
blocked due to a large inter-antidot capacitive coupling. 

Combining the pseudo-spin of a double dot with the 
intrinsic spin, Borda et al. mi predicted an SU(4) Kondo 
effect which has been recently observed [ 51 ] . The use of 
double dots to realize pseudo-spin SU(2) Kondo model 
and its generalizations at v = 1 integer quantum Hall 
regime was proposed in [32] . Here, we extend those ideas 
by studying the realization and transport properties of a 
Kondo impurity coupled to chiral Luttinger liquid edge 
states in the FQH regime. A similar model arises in the 
study of a double dot inserted in a spinless non-chiral 
Luttinger liquid, once the occupancy of the dots is lim¬ 
ited so that they act as an effective spin. In contrast to 
most previous studies that focus on zero temperature, we 
provide analytical expressions for the conductance in all 
asymptotic temperature regimes. In this paper, we only 


deal with fully polarized (or spinless) systems and spin 
refers to the orbital pseduo-spin. 
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FIG. 1: (color online) The system considered here. Two 
antidots enable tunnelling of quasiparticles between the 
outer edge states in the v — \jq Laughlin FQH liquids. The 
gapped incompressible liquid (blue) plays the role of tunnel 
barrier for the quasiparticles. Only one state per antidot is 
considered. The capacitive coupling between the antidots is 
large enough to keep their total relative occupancy constant. 

The model - We consider the system depicted in Fig. [I] 
in which each antidot is represented by a single fermionic 
quasiparticle level. This is valid for small enough antidot 
radius. In this limit the sytem can be described by the 
following Hamiltonian 

H = H 0 - [t R 'iplp, R (0)dR + t L ^l p L (0)d L + h.c.] 

~\~U{d^ R d R + d^ L d L — 1) — tcC R d L + h.c.), (1) 

where th^R^c are tunnelling amplitudes and U is the 
Coulomb energy. Here d L / R annihilates quasiparticles 
on the upper/lower antidot and 'ipqp^R/L^) annihilates 
right/left-moving quasiparticles on the upper/lower edge 
of the Hall bar, with corresponding Hamiltonian Hq. We 
are interested in a parameter regime T,tL,t R ,tc <C D, 
where D ^ min(£e, U ) A. Here 5e is the anti¬ 

dot level-spacing and A is the bulk energy gap (we set 
k R = 1 throughout the paper). Then U limits the an¬ 
tidots charge configuration to (0,1) and (1,0) sectors. 
Sequential tunnelling is blocked in this large U limit 
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and different methods must be developed to study the 
system. Temporarily ignoring the inter-dot tunnelling, 
tc, we see that to transfer one quasi-particle from up¬ 
per to lower edge, we must start in a state with the 
lower dot occupied, pass through a high energy inter¬ 
mediate state with both dots occupied or empty and end 
up with only the upper dot occupied. Thus it is conve¬ 
nient to identify L and R with pseudo-spin up and down 
respectively. A Schrieffer-Wolff transformation [33] then 
yields a Kondo model with impurity pseudospin opera¬ 
tors S = d^aapdp /2 and quasiparticle pseudospin den¬ 
sity J[x) = ip^ p ^ a (x)a a gip qp ^(x)/2. The Kondo interac¬ 
tion, J±[S x J x (0) + S y Jy(0)] + J z S z J z (0) contains Kondo 
couplings J_[_ = 4:tLtR/U , J z = 2{t 2 L -\-t 2 R )/U + SJ Z , where 
5J Z is an additional positive contribution arising from 
the Coulomb interaction between quasiparticles on the 
antidots and edges |[29j [34]. Inter-antidot tunnelling cor¬ 
responds to a magnetic field term in the Kondo Hamil¬ 
tonian, coupled to the impurity spin only, —tcS x . 

While 'ipqp^L /r (t) are simply free chiral fermion fields 
for the integer Hall state occurring at v = 1; for v — 1/3 
it is very useful to bosonize, ip qP: R/L( x ) oc e ±Z(fR / L ^ in 
terms of chiral bosons (Pr/l( x ), obeying the basic com¬ 
mutation relation \ip r/tA x ) , r/tAv) 1 = Eiiriysigntx — y). 
Then Ho = [v/(2irv)] fZ> dx[(d x cp R ) 2 + (d x cp L ) 2 ], where 
v is the quasiparticle velocity. It is then convenient to 
define commuting right-moving spin and charge bosons, 
c p s ,c { x ) = [Tr( x ) A Tl{~ x )]/V^ since only the spin bo¬ 
son appears in the Kondo interaction. Then we obtain 
J- = J x - ijy oc e an q j z — d x (p s /(2 iyvx/2). The 
renormalized Kondo couplings grow larger as the energy 
scale is reduced [35], becoming large at the crossover scale 
T k . For the v = 1 case, T K = De~ x ^ x where A = pj, 
g is the density of states and we have assumed J± = J z 
(similar behaviour occurs in the anisotropic case). For 
v = 1/3, e ly/ ' 2( ' Ps has renormalization group (RG) scal¬ 
ing dimension 1/3. Thus X±(E) = (D/E) 2j/s X j_, so 
T k oc D/X 3 / 2 . 

Conductance - We are interested in the inter¬ 
edge tunnelling conductance, corresponding to back- 
scattering, defined using the charge current operator 
/ = ivetc(d^ L d R — h.c) = 2 vetcS y . In the linear response 
regime the Kubo formula gives m 

G = lim ^ yy ^ \ ( 2 ) 

in units of ue 2 /h where x'yy is the imaginary part of 
the dynamical impurity spin susceptibility of the Kondo 
model, Xyy(u) = -i dte lut ([S y (t), S y ( 0)]). Every 
transmitted quasiparticle contributing to the transport 
involves a spin-flip process at the impurity, relating con¬ 
ductance to the spin relaxation. 

High temperatures: TkRc «f - We may attempt 
to calculate the susceptibility using perturbation theory, 
but this gives a result which diverges at uj —>• 0 : 

x'yy M = - ^ [l7 X 2 Z + Ai (T/D) 2 "- 2 } (3) 


Here 7 ^ and 7 ^ are dimensionless coefficients of 0(1). 
This surprising infrared divergence is not connected to 
the usual renormalization of the Kondo couplings since 
T Tk • Nonetheless, it suggests that an infinite subset 
of diagrams must be resummed to get a finite conduc¬ 
tance [36] • One way to solve this problem is to phe¬ 
nomenologically describe the impurity spin by the Bloch 
equations HZl ESI 

d t {S a ) = [h(t)x(S)] a -^E — a = x,y,z. (4) 

x a 

Here h(t) = (—tc,h y (t), 0) where h y (t) is an infinitesi¬ 
mal time-dependent ^/-component of the magnetic field, 
introduced to obtain Xyy We expect Eq. Q to hold as 
the equation of motion for the averaged impurity spin 
in a theory where the quasiparticles are formally inte¬ 
grated out. Here, xo ~ — 1/4T is the static suscep¬ 
tibility (S )o = Xoh, in the presence of the static field 
—tc- (S a )(t) rotates around the external magnetic field 
and relaxes towards it within time-scale r a because of its 
coupling to the quasiparticles. Therefore, using the def¬ 
inition Xyy{w) = (S y )^ /h y (<jj) for the imaginary part of 
the susceptibility we obtain 


Xyy(v 0 ) 


UXOTy 

1 + t2 C TyT Z ’ 


(5) 


To obtain the conductance, we need r V:Z . The main 
‘Bloch equation’ assumption, justifiable at T , tc, 

is to neglect the frequency-dependence of these rates, 
thus obtaining them from a large frequency limit of our 
perturbative result using Xo/ T z = lim^-^ uox zz {w) oc 
A 2 l (T/D) 2u ~ 2 and y 0 /r y = lim^-**, ux'y y (u) = Xz>? z + 
7 ±A \{T / D ) 2u ~ 2 . So at high temperatures 


G u = i oc 


t 2 

l c 


t 2 

l c 


T 2 (X\ + X 2 ) ’ G " =1/3 K X 2 ± T 2 / 3 D 4 / 3 ' ^ 


(We show explicitly that this result can be obtained, at 
high T, from a resummation of Feynman diagrams in the 
special case X y = A z = 0 in the Supplementary Mate¬ 
rial.) More correctly, A_l, A^ should be replaced by the 
renormalized quantities at energy scale T, but this is an 
unimportant correction assuming T Tk • 

T, Tk tc - In this regime, the impurity spin be¬ 
comes a classical field pointing in the direction of the in¬ 
stantaneous field, h(t) (see Supplementary Material Sec. 
V) so we may approximate: 


H^H 0 T (l/2)Ax|>7x(0) + (h y (t)/tc)M 0)]. (7) 


This corresponds to a direct tunnelling term between 
edges: H T = (l/4)A_ L [e^W/ tc ^ pZ/ (0)^ P ,fl(0) + h.c.}. 
For v = 1 this is a simple non-interacting tunnelling 
model giving a conductance G oc A]_. [More accu¬ 
rately, A_l should be replaced by the renormalized cou¬ 
pling A ±(tc) but this is again unimportant for Tk tc] 
For the fractional quantum Hall case the behaviour is 
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much different [39] since this direct tunnelling interac¬ 
tion is relevant and A_l(T) = (tc/T) 2 / 3 \±(tc)- There¬ 
fore the conductance starts to grow as G oc T -4 / 3 . It 
starts to level off at T of order Tk eventually saturat¬ 
ing at v, corresponding to perfect transmission through 
the double antidots. The nature of this zero temper¬ 
ature infrared fixed point can be straightforwardly un¬ 
derstood from bosonization. The relevant tunnelling 
term, oc — A^ cos[v / 2^ s (0)], pins p s (0) at 0. To under¬ 
stand the physical implications of this boundary condi¬ 
tion note that while the charge boson remains continuous 
at x = 0 at both high T and low T fixed points, imply- 
ing ^fl(0+) - <p L ( 0“) = <Pr( 0“) - <p L ( 0+), the high T 
and low T boundary conditions on the spin boson im- 
ply <p R ( 0+) + = ±[</?ij(0 - ) + <p L ( 0+)], respec- 

tively. Together, these boundary conditions merely im¬ 
ply continuity of Pr/l at the origin at high T but im¬ 
ply (/^(O 41 ) = ~Pl( O 1 ^) at low T, corresponding to a 
breaking of the system into x < 0 and x > 0 parts, 
perfect transmission through the double dots and per¬ 
fect backscattering. The leading low T reduction of 
the conductance is conveniently calculated by consider¬ 
ing the small horizontal current between the nearly dis¬ 
connected x < 0 and x > 0 parts of the system [inset 
of Fig.[2^a)] . This involves electrons tunnelling through 
vacuum (as opposed to quasiparticle tunnelling through 
incompressible liquid) between x < 0 and x > 0 sides 
and corresponds to a term in the effective Hamiltonian 
oc cos[((/?#(0 + ) — pr(0~))/v\, of RG scaling dimension 
1/v. Thus the horizontal conductance, past the double 
dots is Gh oc T 2 T~ 2 = T 4 . By current conservation, we 
expect the vertical conductance through the antidots to 
behave as G v — otfiT/T k ) 4 , for a dimensionless con¬ 
stant of 0(1), a. The behaviour of the conductance when 
Tk < tc for various temperature ranges and v = 1 and 
v = 1/3, is plotted in Fig.[2](a) . The zero temperature 
conductance, as well as the exponens for T tc agree 
with previous numerical results [39] . 

Strong Kondo coupling fixed point, < T < Tk 
- In this parameter regime, the Kondo coupling con¬ 
stants A^ and A^ renormalize to large values but the 
interdot tunnelling, tc , may be treated as a small per¬ 
turbation. The impurity spin is then screened by the 
quasiparticles and, for v = 1, we may apply Fermi liq¬ 
uid theory. The impurity spin, S x appearing in the in¬ 
terdot tunnelling Hamiltonian can then be represented 
by (v/TK)fi>l P (0)cr x 'ip qp (0), the lowest dimension opera¬ 
tor with the correct SU(2) spin transformation properties 
mm- The factor of v/Tk can be inserted by dimen¬ 
sional analysis, recognizing that Tk is the characteristic 
energy scale, or reduced bandwidth at this fixed point. 
The corresponding Hamiltonian is non-interacting, with 
this tunnelling term being marginal under the renormal¬ 
ization group. This leads to the familiar Shiba formula 
[42] giving G oc (tc/T^) 2 . Similar reasoning may be ap¬ 
plied to the is = 1/3 case but now the effective interaction 
oc fi^ p (0)a x fi qp (0) is relevant, with dimension 1/3. Thus 
calculating the conductance to lowest order in tc gives 


(a) 




FIG. 2: (color online) Inter-edge conductance vs. 
temperature for v — 1 and v =1/3 quantum Hall states at 
(a) tc Tk and (b) tc Tk • The dashed lines are 
interpolations which we expect to be qualitatively correct for 
the crossover regimes. Whereas the v — 1 case exhibits a 
crossover only at larger crossover scale max(tc, T^), the 
v =1/3 case has an additional crossover at Tk and T* for 
the case (a) and (b), respectively. Insets: (a) Schematic of 
the stable infrared fixed point in the case of v = 1/3. The 
leading irrelevant processes correspond to electron tunnelling 
through vacuum, (b) The Kondo fixed point conductance in 
the case of v = 1 has a non-monotonous dependence on 
w = tc /Tk, with a peak of 0(1) at w ~ 1. 


G oc (tc/T k) 2 (Tr/T) 2 ^ 1 ^ oc T 4 / 3 . This diverges at 
low T signalling the breakdown of perturbation theory in 
tc- 

T tc Tk regime- For v = 1 there is no sig¬ 
nificant change in behaviour as T is lowered to zero be¬ 
low tc, with the conductance being approximately con¬ 
stant. On the other hand, for v = 1/3, the growth of 
the interdot tunnelling term under renormalization sig¬ 
nals the crossover to the same fixed point discussed above 
for T,Tk <C tc, corresponding to perfect transmission 
through the antidots. Renormalized interdot tunnelling 
becomes strong at the scale T* oc t^/ 2 /T^/ 2 and be¬ 
low this scale the conductance should again crossover to 
v — <a(T/T*) 4 behaviour. The behaviour of the conduc¬ 
tance when tc <C Tk for various temperature ranges is 
plotted in Fig.[2][b). 

Note that the conductance versus temperature looks 
rather similar in the two cases Tk <C tc and tc <C Tk- 
One essential difference is the crossover temperature 
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scales. For v = 1, there is only one crossover which 
occurs at the larger of tc and Tk- For v = 1/3 there 
are 2 crossover scales: and Tk for tc Tk, but 

Tk and t^ 2 /T^ 2 for tc <C Tk- It is also interesting to 
note that, for v = 1 and A^ = A^, the T = 0 conduc¬ 
tance is oc A_l(£c) 2 = 1/ ln 2 (tc/T k) for Tk tc but 
oc (tc/Fx) 2 for tc <C Tk- G(0) decreases as tc/TK be¬ 
comes large or small, going through a peak of 0(1 ) at 
Tk of order tc [see inset of Fig.[2jb)]. 

Conclusion - We have mapped the conductance 
through two antidots in v = 1 integer and v = 1/3 frac¬ 
tional quantum Hall systems onto the susceptibility of a 
Kondo impurity in a Luttinger liquid, analyzed the fixed 
points and calculated the conductance in all asymptotic 
regimes. Calculations of noise, and extension to more 
exotic filling factors v = 5/2 and 12/5 with non-abelian 
statistics are left as future extensions of these results. 
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SUPPLEMENTARY MATERIAL 


In this supplementary material we provide details 
and proofs of some results presented in the paper. Ap¬ 
pendix^ discusses the weak coupling and strong Kondo 
coupling fixed points and the corresponding glueing con¬ 
ditions. Appendix|B| contains some discussion about the 
antidots and experimental considerations. Appendix|C| 
provides a calculation of the susceptibility using the semi- 
classical Bloch equation. Appendix|D] provides the result 
of perturbation theory to second order in Kondo coupling 
but exact in tc- In Appendix|E|we provide a detailed dis¬ 
cussion of the effective Hamiltonian derivation in the case 
of T,Tk <C tc- Finally, Appendix[F| contains an exact 
solution of the case \ y = A 2 = 0, using techniques devel¬ 
oped to study the X-ray edge singularity. We show that it 
is possible to derive the Bloch equation as the high tem¬ 
perature (T >> tc) formula for the susceptibility. The 
exact result in this special case also demonstrates the 
breakdown of perturbation theory at high temperatures. 


for x > 0, in terms of which the Hamiltonian becomes 


H = H 0 + T( 5 '+ e # s (°) + h.c.) 

+Ns z d x e s (o)-t c s x + h z s z (Ai) 

47 TV 

where we also allowed for a detuning h z = 
£r — £l between the two antidots. Here Hq = 
/o°° i^K^s) 2 + (d x 0s) 2 + (d x <t>c) 2 + {d x 9 c ) 2 \ and the 
non-chiral bosons obey [cj) s (x),0 s (y)\ = 2'kiuQ(x — y) 
where Q(x) is the Heaviside step function. 


2. Weak coupling 


The only boundary condition at the weak cou¬ 
pling fixed point is that ip SjC (x) are continuous at 
x = 0, i.e., 6> s , c (0) = 0. Using the definition 

[/ dxp R / L (x),ip qPjR/L (y)] = -ipq P ,R/ L (y) we obtain the 
density operator for both bosons to be Pr/l( x ) = 
2 ^d x (fR/L- From this, assuming a finite length L with 
periodic boundary condition, the mode expansions are 


2ttr 

Tr/l( x ) — — N r /l x + Tr/l, o 


n =1 




±ik n x | -f p 

R/L,n ^ a R/L,n 


Tik n x\-^ 


)e—8-,(A2) 


Here k n = 27 xn/L, and the harmonic bosons and the 
zero mode obey the standard commutation relation 
\- a R/L,n’ a R/L,rJ = 5 nm and [<p R / L ,o,N R/L ] = i. This 
leads to the mode expansion of the charge/spin bosons 


Tc/ s( x ) 


2ixv 


(Nr ± N L )x + 


1 

71 


(Tro T Tlo) + • * *(A3) 


Defining total bulk charge Q = Nr + Nl and spin 2 s z = 
Nr — Nr, for Q even, s z has to be integer, while for Q 
odd it must be half-integer. So, the glueing condition for 
(Q,s z ) at weak coupling fixed point is [41] 


(Q,s z ) = (even,integer) 0 (odd,half-integer). (A4) 


3. Kondo fixed point 


Appendix A: Gluing conditions at fixed points 

1. Folding transformation 

For a discussion of the fixed points, it is convenient 
to fold the chiral bosons, according to 

<Ps c(x) + <Ps c (—x) <Ps,c( X ) - <Ps,c(~ x ) 

K< = - ^^ - 


By power counting, the J z term is marginal but the 
J± terms is relevant. In order to account for this, we 
define dimensionless couplings A^ = 27 tJ z /v and A^ = 
JrD"- 1 and will frequently switch between J Zj r and 
A^jl notations in the following. These couplings grow 
as the bandwidth is reduced [35, and the system flows to 
the Kondo fixed point. Although J z is naively marginal, 
because of the coupling to S z , it controls the scaling di¬ 
mension of the Jr term. This can be seen if we ap¬ 
ply a unitary transformation jMHU H ->■ V/HVp with 
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= exp[i/jLS z (f> 8 {0)] which gives 

H -> Ho + ^(5'+e i( 1 “' i)< M° ) + h.c.) 

+^-^s z d x e s (o) + h zS z 

—tc [S x cos /ufis (0) + S y sin /i(j) s (0)] (A5) 

and changes the dimension of J± from 1 — z/tol — z^(l — 
/a) 2 . In order to understand the strong coupling Kondo 
fixed point, it is convenient to either tune this dimension 
to zero (the so-called Toulouse [45] point) so that it could 
be refermionized [44 or to 1 (the so-called decoupling 
poing) so that it becomes a boundary magnetic field. We 
use the latter approach, which happens at /i = 1 and the 
transverse Kondo coupling becomes J±S X . In the case 
of tc = h z = 0, and if = Ekv, the Kondo coupling 
reduces to a Zeeman field on the isolated (but dressed) 
impurity spin which projects to the ground state of S x 
at low energies. The Kondo temperature is ^ J± in this 
highly anisotropic Kondo model [46]. It is easy to check 
that (p s (x) develops a discontinuity at x = 0 

(p s {x) = V^s{x)V lx = ip s (x) - wis/j,S z sign(x) (A6) 

or equivalently, in the folded basis, the new boundary 
condition corresponds to 0 S (O) = 0 S (O) and a new pin¬ 
ning of 0 S ( 0) = 07tza The pinning is dynamically switch¬ 
ing between these two values as in the instanton-gas rep¬ 
resentation of the Kondo problem [47 . The charge boson 
is unchanged by the unitary transformation followed by 
the projection and the new gluing condition at the strong 
Kondo coupling fixed point [4T is 

(Q,s z ) = {even,half-integer) 0 {odd,integer). (A7) 

This change in the glueing condition, implies that a spin 
boson has decoupled from the edge states to screen the 
impurity spin. We have used the decoupling point to 
discuss the Kondo fixed point, and this requires tuning 
J z to the large value of Attvu which is not physical, as 
the bare Kondo coupling is usually assumed to be small. 
However, it is expected that other values of J z would 
have a similar qualitative behaviour. 


magnetic length £r = but we can assume that 

the antidots are large enough so that the radius of |n) 0 
and \n + 1) 0 are effectively the same [29. and assume that 
charge and neutral sectors decouple. This corresponds to 
the constant interaction model in quantum dots [48]. For 
the antidots, we are interested in a regime where bosonic 
excitation energy Se = 2 ttv/L is much larger than ksT 
and we can assume that harmonic part of the field is in 
its ground state |0) n = 0. Note that |n) 0 and \cp) n 
are analogous to the number of charges and the excited 
states of normal quantum dots. 


We also need to take into account the Aharonov- 
Bohm contribution of the magnetic flux going through 
the antidots. The number of quasiholes Ur in antidot R 
is such that it is equal to the number of flux quanta going 
through Pr/P o = ur where <po = h/e. This is another 
way of stating that Rr = ^/2ur£b, at the ground state. 
These numbers change as we change pR. This is done by 

by E c (n R - 4>R/<t> o) 2 with 


replacing Ecn 2 R in Eq. (B1 


ipo = h/e. We are interestec 


in a regime where two states 


in the same antidot with ur and ur + 1 quasiparticles 
become degenerate. This is possible for Pr/P o = 2mR+l 
where tur is an integer. Also we need similar degeneracy 
to be valid for the second antidot <Pl/4>o = 2 tyil + 1. To 
have both of these at the same magnetic field, we obvi¬ 
ously need some fine tuning of the area of at least one of 
the antidots. We assume that this is possible by tuning 
the voltage applied to the gates that defined the antidots 
at the first place, or by a combination of voltages applied 
to the outer edges. To capture deviation from this ideal 
case one can add a term h z S z to the Hamiltonian, but 
we assumed a perfect tuning in the paper. 

If the temperature is low enough (T Ep) 

so that the bosonic modes of the antidots are not excited, 
they effectively behave as hardcore fermions [29.. To see 
this, following m we assume that N and N +1 states of 
the antidot are degenerate and denoting them by |0) and 
11), it can be seen that due to the commutation relation 
[upo ,7V] = i, the operators s± oc e ±lcpo are raising and 
lowering operators of the ’’spin” made of |0), |1). From 
[ifo, A] = i, it follows 


[TV, e± i(po } = ±e ±icpo and Ne ±iipo = e ±iip0 (TV ± 1) 


Appendix B: Antidots 


We can find the s pect rum of the antidot by inserting 
the mode expansion (A2) into the free Hamiltonian 


H = 


2nv 


f L 

/ dx{d x cp) 2 = vka\.a k 0 E c n 2 , (Bl) 
k> o 


where Ec = 27 xvv/L acts like the “charging energy” of 
the antidot and L is its circumference. We see that the 
number of quasiparticles n is a good quantum number 
ft \n) Q = n\n) Q . These two sectors are coupled to each 
other because L = 2ttR and R = %/2N£r in terms of 


These can be combined with the bosonization Klein fac¬ 
tors T l,r to represent the creation and annihilation op¬ 
erators for the additional fermion on the dot. 

4 = T L e iVL0 and 4 = T ( B2 ) 


Appendix C: Bloch equation - Non-zero tc 


Considering that h = {—tc, h y , 0) and X y ^ X z , there 
is no spin symmetry present and we have to allow for 
different relaxation rates along each direction. Therefore, 
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we can write the Bloch equations [Eq.(4) of the paper] as 
dt(S x ) = h v (S z )- (Cl) 

T~X 

dt < Sy ) = tc (s z ) - (C2) 

T y 

dt (S z ) = -tc (Sy) - hy (S x ) - (03) 

T z 

where (S' 0 ) are the components of the steady state mag¬ 
netization. To find the steady state magnetizations, we 
do a rotation (tail'd = — h y /tc ) 


we expect the second term to be 0(h 2 ). Then everything 
simplifies: We get (S x ) u ~ —27r5(u))xotc, he. constant in 
time. Thus the convolution in the last line also simplifies 
and we obtain 


(1 - iur z ) (S z ) u « -t C T z [ ( Sy ) w - Xohy(u)\ , (CIO) 
from which we get 

v (,,)= lj m Xo(l -iUT z +t%T z Ty) 

yy hy^O hy{uS) (1 — iuJT z )(l — lUJTy) + tgTzTy 

(Cll) 


with the imaginary part 


OW=(“) (C4) 

v ' \ smd cos'd J v 1 

on 

(h x 0) = (-tc h y ) 0(t9), (|“) =0(^)(|) 

to obtain h x = cos $(—tc + hy/tc) = —+ 0(h 2 ). The 
Hamiltonian is diagonal in this ‘tilde’ basis and we find 

(S' 0 ) = cos'd ^ tanh —/A cos'd, (C5) 

= _ ^ tanh ^ sin ^’ (>s°) = o ( C6 ) 

Since eventually we are interested in Xyy = 

d (S) y /dh y \h y =o we can drop 0{h 2 ) and the above 
results simplify to 

(S^itaah^, (C7) 

(SJ>=0 (CS) 

For T tc where we expect the Bloch equation ap¬ 
proach to be valid, tanh(tc/d/2) tttc/3/2 and we get the 
linear response result (S°) w xo^y 5 but for T <C tc, we 
have tanh(tc/d/2) ~ 1. More generally, we can define an 
effective static susceptibility given by (S' 0 ) ~ xoh , where 

Xo (T) = —— tanh = xo (t —)>- ^ ) (C9) 

V ' 2t c 2 AU V 2 tanh ^ V ' 

Fourier transforming, we obtain 
(1 - iu>T x ) (S x ) u = -27T<5(w)xo ff fe + T x h y (u) * (S z ) w 
(1 - iuTy) ( S y ) u = xfh y (uj) + t C T y (S z ) u 
(1 - t z ) (S z ) u = -t c T£ (S y ) u - T z h v (u) * (S x ) u 

These are easily generalized to the memory-full case, by 
allowing a frequency-dependence for T a (u). These set of 
equations are difficult to solve. One approximation that 
greatly simplifies this, is to ignore the fluctuations of the 
spin along the external field, x direction. This amounts 
to dropping the second (convolution) term on the right 
hand side of the first equation and makes sense because 


,// / a = ‘'-W# + t%T z T v + uCf\ 

VV [^ + (t 2 C -U 2 )T z T y \ 2 +U 2 (T y + T z ) 2 


(C12) 


Let us look at this formula, in various limits. Without 
Kondo coupling T x ^ y ^ z —>> oo, and we get 


Xyy(u) - 


tcX o 

“ 2 ~tc 


tc tanh(ffic/2) 
2 uj 2 — t 2 c 


(C13) 


For t Cy jT z Ty < 1, we basically get the simple result 
X yy (w) = Xo/(l — iuTy) that we would get if we had 
neglected tc from the beginning. Generally, we see that 


Xyy 0 ) 


UXQTy 
1 + t 2 C TyT Z ’ 


But for uj oo 


-8^cX0 Ty 

1 + t2 C T V T z 


(C14) 


XyyC -> OO) = — = lim U>Xyy(uj) (C15) 

yy UOTy Ty OJ^OO yy 

To obtain the conductance, we also need r z which can be 
obtained using 


— = lim ux zz (u) (C16) 

T z w—>-oo 

or from r y with a rotation along S x , i.e. by interchanging 

Xy ^ X Z . 


Appendix D: Susceptibility to order 0( A 2 ) but exact 

in tc 

In this section we provide the result of perturbative 
calculations of the imaginary part of the susceptibility to 
second order in Kondo coupling but exact in inter antidot 
tunneling tc • The goal of this section is to demonstrate 
that once a finite tc is included the infrared divergence 
of the perturbation theory is cut off. Using equation of 
motion techniques, it can be shown that the correlation 
functions to second order in Kondo coupling are 

*T = %iX5sr lm[n > )] m 
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Here Hf» (, J a Ja)u are retarded correlation func¬ 

tions of the current operators, 


u^ y (u,p = i) = n *(u,u = i), n£(w,i/) = 

n %y(u,v < 1/2) = -(y) sin(7T v)b(v- ^,1-2^ 


where B(x,y) is the beta function. The X 2 contribution 
has a more complicated form 

Xyy(u < tc) = ~g{^ ~ 2 + 


--tanh(/3ic/2) 


where he function iJ„(w) for v = 1 is 


H u= \{uj) = Im 



iWgtal) ^( 1 + 


i/3(t c ~uj) \ - 
2tt y 


CJ + tc 


UJ ~t C 


in terms of digamma function 'ip(z). Note that in the 
limit of tc —> 0, these results reduce to Eq. ^ of the 
paper. 


Appendix E: Effective Hamiltonian for T, Tk tc 

Temporarily igoring the infinitesimal time-dependent 
part of the field, h y (t ), introduced to calculate the dy¬ 
namical susceptibility, y" , it is clear that in this regime 
we may replace S x by 1/2 (and S y,z by zero) since the 
impurity spin is polarized by the strong field. Now con¬ 
sider the effect of h y (t). We again wish to integrate out 
the impurity spin to obtain an effective Hamiltonian for 
the quasiparticles. It is now not appropriate to consider 
ant relaxation terms in the Bloch equations, since such 
terms arise from the integrating out the quasiparticles in¬ 
stead. So, let’s consider the solutions of the simple spin 
torque equation 


dtS = h(t)xS(t) (El) 

with h(t) = (— tc , h y (t), 0), h y (t) = etc eoscc/t, taking the 
limit where both e 0 and u 0. We also assume that 
the oscillating component of the field is turned on slowly 
in the infinite past. Thus we write: 

S(t) = (1/2) (1,0,0) — S'. (E2) 

We will see that S' is O(e). Working to first order in e, 

h x S/t c ~ — (1/2) (0,0, e cos ut) + (0, -S z ', S y ') (E3) 


Thus 


1 dS z ' 
tc dt 
1 dS y ' 


S yf — - cos ut 
2 


-S z '. 


(E4) 

(E5) 


Thus 


1 d 2 S z ' 
tc dt 2 


—tcS zf + -ecu smut. 


(E6) 


with solution 
S z ' 

Thus 


with solution 


et c uj 

2(^-^ 2 ) 


sin ut 


eu . 

-— smut 
2 tc 


dS y ' eu . 

—;— ~ — —si 
dt 2 


S y ' = ^ COS Ut. 


(E7) 


(E8) 


(E9) 


In summary 

S ~ (1/2)(1, —ecosut, —e{u/tc) smut). (E10) 
At small u we may drop the last component, giving 

S^-(l/2t c )h{t). (Ell) 


Thus we see that, even with no relaxation term in the 
Bloch equations, purely a precession term, the spin 
tracks the instantaneous time-dependent field, in the 
limit where the time-dependent term in the field is small 
and slowly varying. 


Appendix F: Exact spin susceptibility when 

A y A z — 0 

1. Derivation of Bloch equation result for tc ~ 0 

The fact that Eqs. (3) and (6) of the paper contain 
a simple summation of A^ and X 2 contributions to low¬ 
est order, suggests that the same IR divergence and the 
necessity to use Bloch equation occurs even if only one 
Kondo coupling, say A^ ^ 0, is non-zero. In that limit 
and assuming we could neglect tc ~ 0 at T tc, it is 
possible to map our Kondo problem to the X-ray absorp¬ 
tion problem and find a non-perturbative formula for the 
susceptibility. 

With only A 2 nonzero, the Hamiltonian is H — Hq + 
X z S z J z (ti) and the Kondo interaction is just a boundary 
magnetic field, depending on the spin state of the im¬ 
purity. We are interested in the dynamic susceptibility 
(. S y S y ) UJ1 defined after Eq. (2) of the paper. But S y it¬ 
self, is not present in the Hamiltonian and its influence 
is just suddenly switching the sign of the boundary mag¬ 
netic field, via S y |t) = i ||) and S y ||) = — i |t)- Before 
this switching, the spin-up fermions see a phase shift and 
spin-down fermions another and these two phase shifts 
are suddenly switched. The ground states before and 
after switching are orthogonal to each other in the ther¬ 
modynamic limit m and the transition creates lots of 


tc dt 

















electron-hole pairs, the so-called orthogonality catastro¬ 
phe. It is more convenient to discuss this in terms of the 
spin-up/down bosons. Setting the velocity to 1, dropping 
constants and following [12 by introducing = Pr/l ? 
the Hamiltonian density at v = 1 is 

n = n 0 + 6(x)^s z (d x ^-d xn ), (fi) 

where Ho = (d x (p^) 2 + (d x ip±) 2 . Depending on the state 
of the impurity spin |t) or ||), the Hamiltonian breaks 
into two sectors H = H+ |t) (t| + H_ \\) (||, where using 
= gJ z = J z / we get 

n± = [d x ^ ± ^ 5{x .)] 2 + [0 xn T ^(x)} 2 , (F2) 

up to a constant. Defining A' = A^/4, these two Hamil¬ 
tonians are related to Ho by the Schotte-Schotte unitary 
transformation l52l 


that 


U = e -^Vt(°) e + iA Vr(o) (pg) 


H+ = U^HoU, H- = UU 0 U\ (F4) 


e 4i7rA ' 2 ( e 


-2iA'[y t (i)-y t (0)]^ ^ e 2jA'[v3j.(t)-V4.(0)]^ 


0 4i7ruX'^ 


^8uX' 2 


Doing similar procedure for the other terms and summing 
up all the terms, at zero temperature we obtain 




/ A ©ft) . 

(*) = -^ sm7T 9, 


g = 4A' 2 = Y (F 8 ) 


This correlation function is a power-law in absence of 
any bulk interaction, because of the physics of orthogo¬ 
nality catastrophe [EEJE 2 - Using conformal mapping to 
a finite-radii cylinder t ^ sinh ^, we can bring this 
retarded function to the finite temperature, 


X 




tt\ 2 9 Q(t) sin7 ig 


sinh 1 j 


2 g m 


(F9) 


The Fourier transform 

[2ll\ 2 9~ l 


gives 


Xzz(v) = ~{-p ) B(g -iui/3/2w,l-2g) 


sin ng 


We are interested in the limit of small Kondo coupling, 
g 0. Therefore, using properties of the beta function 


Writing S y as 2iS y = S + — S and applying the unitary 
evolution operator 

2 iS y (t) = e itH +S + e~ itH - - e itH ~S~e~ itH +. (F5) 

Here, we used that before/after applying S + , the system 
has to be in the —/+ sectors, respectively. Inserting this 
into dynamic susceptibility and dropping S ± S ± terms, 
we obtain 

X R (t) oc 0(t) [ ( e itH +e~ itH -S + S~) + (e itH ~e~ itH + S~S + ) 
- (e itH +e~ itH -S~S + ) - ( e itH -e~ itH +S + S~) 


R 

A,zz 


B{g-ioj/3/ 27r,l)-y = 


Xo 


1 — IUJTk 


Using Eqs. (F4) we can write e ±ltH + = U^e ±ltH °U and 
e ±itH- _ jJ e ±itH o jj f . \y e p ave apply the same pro¬ 
cedure as in Eq. ( |F5| ) to the Boltzman factors. If the 
correlation function contains the Boltzman fac¬ 

tor becomes e~^ H jZ e~P H + /Z = e~^ H °U/2Zo and 
similar version for the S~ S+ terms. After these substi¬ 
tutions, the spins have done their job and can be simply 
dropped from the correlation functions and we arrive at 

X zz (t) oc G(t) (e iHot U 2 e- iHot U^) + {e iHot U ]2 e~ iHQt U 2 ) 

- (U^ 2 e iHot U 2 e- iHot ) - (U 2 e iHot U t2 e- iHot ) 

oc 0(t) [ (U 2 (t)U t2 ) + (U t2 (t)U 2 ) 

- {U ]2 U 2 (t)) - {U 2 U ]2 {t)) ]. (F6) 

The first term (U 2 (t)U^ 2 ) at T = 0 reduces to 

e -2iAVtd) e 2iA Vt(°)^ ^ e 2 * A Vrd) e -2« A Vr(o)^ 


where yo = —&/4 and P/tk = 27 ig = 2n\ 2 /A. We have 
dropped g in the second argument of the beta function, 
but kept it on the first argumen t. Th is is what one would 
obtain assuming tc ~ 0 in Eq. (Cl 1) which relies on the 
Bloch equation approach. This form was proved using 
some analytical assumption [37| or by using the phe¬ 
nomenological Bloch equation [38] . but we provided an 
exact derivation here, when \ x y = 0. We expect similar 
result for A^ ^ 0 but X Zjy = 0 by a spin rotation along y 
direction. 


2. Non-zero tc 

When tc is non-zero, for \ y ^ 0 or A 2 ^ 0 the Hamil¬ 
tonian contains non-commuting spin terms and the prob¬ 
lem is complicated. However, the special case of only 
A x 7 ^ 0 (but \ VrZ = 0) can be still solved exactly using 
techniques similar to those described above. Note that 
having A^ ^ \ y is unphysical as a Schrieffer-Wolff trans¬ 
formation would always produce equal transverse Kondo 
couplings. Nevertheless, this unphysical case can be used 
as a check on our Bloch equation result. Following the 
same technique as in previous section, it can be easily 
shown that 


( t) = __r 

Xyy{> 4[£ sinh(7rf//3)p l 


e ~ itct A - e itct A 


where 


A = — cos(Tig) tanh(7r6) + i sin(7r g). 


*}, (F10) 


(FH) 


(F7) and b = /3tc/2n. The Fourier transform of (F10) gives 
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FIG. 3: (color online) The special unphysical case of 
X y = A z = 0, but Acc ^ 0 and tc ^ 0. Conductance G(b,g) as 
a function of T/tc for various values of g = A^ /4 on log 
scale. These values are g — 0.002, 0.004, 0.006, 0.008, 0.010, 
and 0.012, corresponding the lowest to highest conductance. 
The exact result is compared with perturbation theory to 
second order in Kondo coupling and exact in tc (red color) 
and the Bloch equation result (green). 


-4ix£,M = Gn 


/ 27 T \ 2 £ _1 ' 


\pj 

-a*b(. 


AB 


(•- 


i/3(u> - t c ) 


27r 


9 


iP(u + t c ) 

2tt 


, 1 - 


,1-25) 

25)](F12) 


Using Eq. (2) of the paper, the conductance can be writ¬ 
ten as a closed formula 


G = —47r6 2 Im j A(5, g)B(g + ib, 1 — 2 g) 

x [ip(g + ib) -ip(l-g + ib )] |(F13) 


where again 'ip(z) is the digamma funciton. Note that the 
conductance is a function of b = fitc/^ and g = A^/4 
only. Higher values of Kondo coupling squared g corre¬ 
spond to larger conductance. This function is plotted in 
Fig.[3j as a function of {2nb)~ 1 = T /tc for various values 
of g, and it is compared with perturbation theory result 
[Eq. @] and the Bloch equation result [Eq. ( |C14 )]. Al¬ 
though the second order perturbation theory (but exact 
in tc) is sufficient at T <C tc, it fails in the opposite 
regime of T tc as pointed out in the paper. On the 
other hand, the Bloch equation result provides an accu¬ 
rate estimation of the conductance in this high tempera¬ 
ture regime of T >> te¬ 
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